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Abstract. We use the transverse-momentum dependence of the cross section for diffractive dissociation of 
high energy pions to two jets to study some non-perturbative Light-Cone wave functions of the pion. We 
compare the predictions for this distribution by Gaussian and Coulomb wave functions as well as the wave 
function derived from solution of the Light-Cone Hamiltonian in the Singlet Model. We conclude that this 
experimentally measured information provides a powerful tool for these studies. 
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1 Light-Cone Wave Functions 

One of the most interesting subjects of particle and nu- 
clear physics is understanding the internal structure of 
hadrons. It bears directly on the fundamental interac- 
tions of quarks and gluons that create the hadronic bound 
state. It is also an essential ingredient in understanding the 
hadronic strong, electromagnetic and weak interactions. 
A very powerful description of the hadronic structure is 
obtained through the light-cone wave functions. These 
are frame-independent and process- independent quantum- 
mechanical descriptions at the amplitude level. They en- 
code all possible quark and gluon momentum, helicity 
and flavor correlations in the hadron. The light-cone wave 
functions are constructed from the QCD light-cone Hamil- 
tonian [1]: H2c D = P + P~ - PL where P ± = P° ± P z . 
The wave function iph for a hadron h with mass sat- 



isfies the relation: ,\i/> h ) = M£\ip h ). 

The light-cone wave functions are expanded in terms of 
a complete basis of Fock states having increasing complex- 
ity. In this way the hadron presents itself as an ensemble of 
coherent states containing various numbers of quark and 
gluon quanta [2]. For example, the negative pion has the 
Fock expansion: 



> k -L* ) I ud ) + /7 >i,ku)| udg) 
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representing expansion of the exact QCD eigenstate at 
scale A in terms of non-interacting quarks and gluons. 
They have longitudinal light-cone momentum fractions: 



p+ 



u0 



and relative transverse momenta 



i 



= 0, 



(1) 



(2) 



The form of ip n / H (ui,hxi) is invariant under longitudi- 
nal and transverse boosts; i.e., the light-cone wave func- 
tions expressed in the relative coordinates Uj and fc^i 
are independent of the total momentum (p + , pj_) of the 
hadron. The Fock states are off mass shell with masses of 
M n = Y^i=i{k\i + m 1)l u i where rrii are the quark (cur- 
rent) masses. The first term in the expansion is referred to 
as the valence Fock state, as it relates to the hadronic de- 
scription in the constituent quark model. The higher terms 
are related to the sea components of the hadronic struc- 
ture. It has been shown that once the valence Fock state is 
determined it is possible to build the rest of the light cone 
wave function [3,4]. This was done for the pion using Dis- 
cretized Light-Cone Quantization (DLCQ) on transverse 
lattice [5,6]. 

The hadronic distribution amplitude <p(u, Q 2 ) is the 
probability amplitude to find a quark and antiquark of 
the respective lowest-order Fock-state carrying fractional 
momenta u and 1 — u [7,8]. The pion distribution ampli- 
tude and the light-cone wave function of the respective 
Fock state ip are related through [7, 8] : 



(3) 
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Fig. 1. Fits to the Singlet-Model wave function, (a) is the fit 
to the full range, (b) is the fit to the low range first term, and 
(c) is the fit to the higher range, second term. The resulting 
parameters are: a = 0.7, p a = 0.515 GeV/c, b = 2.55, pb = 
1.58 GeV/c. 



Q 2 



u(l — u) 



(4) 



Two functions have been proposed to describe the momen- 
tum distribution amplitude for the quark and antiquark in 
the \qq) configuration. The asymptotic function was cal- 
culated using perturbative QCD (pQCD) methods [7-11], 
and is the solution to the pQCD evolution equation for 
very large Q 2 (Q 2 — > oo): 



4>as(u) = VSu(l - U) . 



(5) 



Using QCD sum rules, Chernyak and Zhitnitsky [12] pro- 
posed a function (CZ) that is expected to be correct for 
low Q 2 : 

<j> cz {u) = 5\/3u(l -u)(l - 2u) 2 . (6) 

In recent experimental work [13] it was concluded that the 
asymptotic distribution amplitude describes the pion well 
for Q 2 > 10 (GeV/c) 2 , see also Fig. 2b. 

Alternative approaches are based directly on the wave 
function, particularly for the non-perturbative low Q 2 re- 
gion. Jakob and Kroll [14] proposed a Gaussian function: 
tp ~ e~P k± . Pauli [15] considers the L z = S z = com- 
ponent of the ud wave function in the Singlet Model: 
iP(uM) = ^„j(«,kL;U). Then H$j D \i/; h ) = M 2 |^> 
translates to Eq.(8) of [15]: 



TTt ~\~ k | OL 

— ±- ^(«,k ± ) - —j x 

u(\ — u) 3tt z 



M 2 ^(w,kj_) = 

r du'd 2 ^ ^K,kj) 

X J Ju(l - u)u'(l - u') \Q 2 V 2 + Q 2 



(7) 



Am 2 



2fi 2 



for equal masses mi = TO2 = m, and for the mean Feyn- 
man four-momentum transfer Q 2 of the quarks. Replacing 



the integration variable u by k z according to 



u = 



1 + 



v/m 2 + k 2 + fc 2 



Inversely, one expresses fc z by w with 
fc 2 = (m 2 + kj 2 ) ^" ~ ~ 



u(l — u) 



(8) 



(9) 



The substitution allows to introduce the 3-vector p = 
(kj_, k z ) and to trans-scribe this integral equation into 



[4m 2 4 
4m 2 

~Q 2 ~ 



4p 2 ] ^(p) 

% 2 + Q 5 



2a f d 3 p' (p(p') 



3vr 2 7 my/AtpjA{p 
= MV(p) . 



x (10) 



This equation was solved for i^s(p) numerically [15], in the 
non-relativistic approximation A(p) = y 1 + p 2 /m 2 <~ 1. 
In this approximation, the mean four-momentum becomes 
Q 2 = (p p ') 2 , and the light-cone wave function tp(u, kj_) 
is related to the reduced wave function ip(k z ,k±) by 



y(fc z ,kj_) 



(11) 



for details see [15]. The parameters a = 0.69, m = 0.406 
GcV, and /i = 1.33 GeV, yield the correct mass-squared 
eigenvalue M 2 = (140 MeV) 2 for the pion and all other 
iso-scalar mesons. 

The isotropic numerical wave function <p(p) is parame- 
trized [16,17] for p < 0.9 GeV/c as: 



¥>(P) 



i2 ' 



p a = 0.515 GeV/c , (12) 



that is, like a Coulomb wave function with an abnormally 
large Bohr momentum: p a > m. In Fig. 1 we present a 
parametrization of the full range of the calculated wave 
function [17] by a two-term Coulomb wave function. The 
first term is kept as in Eq. 12 with p a — 0.515 GeV/c. 
The mean momentum of the second term, pb, and the two 
normalization constants a, b are free parameters: 



+ 



(13) 



By extending the parametrization to the higher momen- 
tum range we can expect this function to be relevant to 
the measured k±_ distribution. 

The light cone wave function (Eq. 11) with the variable 
transform (Eq. 9) becomes 



ip(u,k.j_) = 2Af 



a[4u(l-u)] 3 b[Au{l-u)Y- 
H + K 2 a {u)f 



H + k 2 {u)Y 

with K 2 (u) = 4w(l - u)p 2 a + (2u - ifm 2 , 

Kl{u) = 4u(l - u)pl + (2u - ifm 2 . (14) 
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Fig. 2. Comparison of the kt dis- 
tribution of acceptance corrected 
data with fits to cross section de- 
pendence (a) according to a power 
law, (b) based on a nonperturba- 
tive Gaussian wave function for low 
kt and a power law, as expected 
from perturbative calculations, for 
high kt. Reproduced from Ref.[13]. 



The main experimental test of the wave function in the 
non-perturbative regime has traditionally been through 
measurements of the pion electromagnetic form factor [18]. 
However, as shown by the authors, there is little sensitiv- 
ity of the form factor to the distribution function. It is 
therefore desirable to develop alternative ways for these 
tests. 



2 Measurement of the Light-Cone Wave 
Function 

The recent measurement of the pion light-cone wave func- 
tion [13] is based on the following concept: a high energy 
pion dissociates diffractively on a heavy nuclear target. 
The first (valence) Fock component dominates at large 
Q 2 ; the other terms are suppressed by powers of l/Q 2 for 
each additional parton, according to counting rules [18, 
19]. This is a coherent process in which the quark and 
antiquark break apart and hadronize into two jets. If in 
this fragmentation process the quark momentum is trans- 
ferred to the jet, measurement of the jet momentum gives 
the quark (and antiquark) momentum. The proportional- 
ity of the differential cross section with respect to the jet 
momentum to the distribution amplitude (squared) was 
asserted [20]: 



<jn 



du dMj ■ d 2 P Nt 



= 2.6 GcV~ 6 



GeV 



<m (15) 



and Umeasured = „. - ■ h is the measured trans- 

verse momentum of each jet. It is assumed that k t (jet) ~ 
k±(quark). This relation was also studied via Monte carlo 
simulations in order to verify the proportionality and take 
into account smearings in the fragmentation process and 
kinematic effects [13]. From simple kinematics and assum- 
ing that the masses of the jets are small compared with 
the mass of the di-jets, the virtuality and mass-squared 

of the di-jets are given by: Q 2 



M 2 DJ 



. The 



1/(1— u) 

results of the measurement show that for k t >1.5 GeV/c, 
which translates to Q 2 > 10(GeV/c) 2 , there is good agree- 
ment between data and the asymptotic wave function. The 



conclusion was that the pQCD approach that led to con- 
struction of the asymptotic wave function is reasonable 
for Q 2 > 10(GeV/c) 2 . At lower values contributions from 
non-perturbative effects may become noticeable. In this 
work we focus our attention to this transition region. 



3 The /^-Distribution 



The kt dependence of diffractive di-jets is another observ- 
able that can show how well the various wave functions 
describe the data. As shown in [20] this dependence is 
expected to be: ^ <~ k^ 6 for one gluon exchange pertur- 
bative calculations. The experimental results are shown in 
Fig. 2, reproduced from Ref. [13]. In Fig. 2(a) the results 
were fitted by fc™ for k t > 1.25 GeV with n = -9.2 ± 0.5. 
This slope is significantly larger than expected. However, 
the region above k t ~ 1.8 GeV/c could be fitted (Fig. 2(b)) 
with n = —6.5 ± 2.0, consistent with the predictions. This 
would support the evaluation of the light-cone wave func- 
tion at large k t as due to one gluon exchange, as is the 
asymptotic wave function. 

The lower fc t -region can be considered as a transition 
from the perturbative to the non-perturbative regimes. 
The experimental results go down to kt ~ 1 GeV/c, still 
large to be considered as a fully non-perturbative regime. 
In [13] this region was fitted with the non-perturbative 
Gaussian function: tp ~ e~P kt [14], resulting in /3 = 1.78± 
0.1. Model-dependent values in the range of 0.9 - 4.0 were 
used [14]. This fit, although resulting in the parameter (3 
being consistent with theoretical expectations, is not very 
satisfactory. As seen in Fig. 2(b) the curved shape of the 
theoretical calculation is not observed in the data. 

We now turn to comparison of the experimental data 
with the prediction of the wave function based on the Sin- 
glet Model (section 1). For this purpose we begin with the 
same expression used by the authors of [13], 



cx \a s {k 2 )G{u,k 2 t )\ 



(16) 
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Fig. 3. Comparison of the kt distri- 
bution of acceptance corrected data 
with fits to cross section derived 
from (a) two-terms wave function 
representing the full range of the 
Singlet-Model wave function, (b) 
Two-term wave function for low kt, 
and a power law, as expected from 
perturbative calculations, for high 
k t . 



with a s (k 2 )G(u, k 2 ) ~ fc t 2 , which is based on factorization 
and which actually is a double differential cross section. 
Using it with the from Eq. 14 gives 



d 2 * =M Ml - ")] 3 
du dkt : ' ' 



1 - 



2 k 2 



1 



km 

k 2 t 



(17) 



up to a constant A/i, and up to similar terms with K 2 . 

The fct-distribution is obtained from Eq. 17 as a single 
differential cross section by integrating over u. But since 
d 2 a/du dk t is strongly peaked at u = |, an exact treat- 
ment is both complicated and unnecessary. Rather, it is 
approximated by removing the slowly varying terms from 
the integral, 



da 

dk t 



1 1 (gg) 
2 k 2 



1 + 



Jo 



[w(l 



(18) 



and similarly for the &-terms. The average value (K 2 ) is 
evaluated in two different approximations: 



<*2> = Pi, 

8 1 

and (Kl) = -pi + -m 2 



(19) 
(20) 



and similarly for the &-terms. The former is obtained by 
the peak value {K 2 ) = K^{\), the latter is obtained by 
the weighted average: 

wi-» = f> f ,"?"°' 1 ,"" = l. < 21 ' 

Jo duw(u) 9 

L du w(u) <2u - l) 2 1 , 
((2u-lf) = Jo \> \ ' =-, (22) 

Jo y 

with the weight function w(u) = [m(1 — u)] 3 . For the par- 
ticular case p a = m both approximations coincide. 



In Fig. 3(a) we compare the prediction of Eq. 19 with 
the data over the whole measured range. Only the total 
normalization constant ATi is a free parameter. The quality 
of the fit (x 2 = 1-7) improves significantly when we repeat 
the fit only in the low k t region (x 2 = 0.8). This is shown 
in Fig. 3(b) with the higher k t range left as in Fig. 2(b). 



4 The ^-Distributions 

Finally, we consider the w-dependence of the diffractive di- 
jets by integrating the double differental cross section (Eq. 
17) over the transverse momentum k t . E791 [13] measures 
conditional u-distributions 



da 
du 



with ki < kt < k u 



(23) 



(ki,k u ) 



for the momentum intervals 1.25 GeV/c < fc* < 1.5 GeV/c 
and 1.5 GeV/c < k t < 2.5 GeV/c. For comparing that 
with theory one should integrate Eq. 17, and we do it in 
terms of the auxiliary function G(u) = G(u; ki): 



G(u;fc,) 



d 2 k t 



du dk t 
7rA/i [u(1 — u 



Qikl - kf) 



(24) 



)f jdz 



[z+ Kl{u)\ 



The integration can be carried out analytically: 

G( U ;fc,)=7rM ^-^\ x 
[k 2 + Kl{u)] 7 



(25) 



— - ^K 2 a (u) + ^K 4 a (u) + —Kt{u) 
4 20 aW 40 aW 280 aV ' 



The conditional u-distributions become then 
= G(u;k t ) - G(u;k u ) . 



da 
du 



(26) 



(fcl.fcu) 

In Fig. 4 we compare the predictions of Eq. 26 to the 
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Fig. 4. The asymptotic distribution [u(l — u)] 2 (solid line) is 
compared with two u-distributions calculated for 1.25 GeV/c < 
k t < 1.5GeV/c, and 1.5GeV/c < k t < 2.5GeV/c which 
almost coincide (dashed lines). All curves are normalized at 
tt=0.5. 

asymptotic distribution [u(l— u)] 2 of [7-11]. Only the over- 
all normalization is used as free parameter. They are nar- 
rower than [u(l — u)] 2 mostly due to the factor [u(l — u)] 3 
in Eq. 24. 



5 Discussion and conclusions 

The agreement of the calculated ^-distribution with the 
data in the transition region of Fig. 3 is interesting since 
the wave function is a non-perturbative solution of a rela- 
tively simple model for the Light-Cone Hamiltonian. The 
model was developed for the bound state problem in phys- 
ical mesons [15], where the low k± -properties are impor- 
tant. It is surprising that this model describes also the 
large /c^-properties in the tail of the wave function, with- 
out that a re-fit of the parameters is necessary. 

The non-relativistic approximation to Eq. 10 predicts 
a u-distribution quite similar to that predicted by the 
Asymptotic wave function. This is a result of it being 
a solution of the Light-Cone Hamiltonian for a \qq > 
system that conserves angular momentum. Other non- 
perturbative wave functions such as that derived from a 
harmonic oscillator potential do not have this property. 
The agreement of the predicted u-distrubution with the 
E791 data [13] is not as good as that of the Asymptotic 
wave function but this is not surprising as these data were 
taken in a kt region that is in the k± tail of this function. 

The non-relativistic approximation to Eq. 10 is a tech- 
nical simplification but not a compelling part of the model. 



It can be relaxed in future work, without major difficul- 
ties. 

In general, we can note that the fcj_-distribution is a 
powerful tool for studying wave functions in the pertur- 
bative regime and in the transition region between the 
perturbative and non-perturbative regimes. 

The data of E791 were taken more than ten years ago. 
With realistic light-cone wave functions now existing or 
coming up in the forseable future, dedicated diffraction ex- 
periments, possibly extended to kaons and other hadrons, 
should be very valuable. 



References 

1. S.J. Brodsky, D.S. Hwang, B.Q. Ma and I. Schmidt 
Nucl.Phys. B593, 311 (2001). 

2. S.J. Brodsky, SLAC-PUB-7870; hep-ph/9807212. 

3. A.H. Mueller, Nucl. Phys. B415 (1994) 373. 

4. F. Antonuccio, S.J. Brodsky, and S. Dalley, Phys. Lett. 
B412 (1997) 104. 

5. S.J. Brodsky, H.C. Pauli and S.S. Pinsky, Phys. Rep. 301, 
299 (1998). 

6. S. Dalley, Nucl. Phys. B (Proc. Suppl.) 90 (2000) 227. 

7. G.P Lepage and S.J. Brodsky, Phys. Rev. D22 (1980) 
2157. 

8. S.J. Brodsky and G.P. Lepage, Phys. Scripta 23 (1981) 
945. 

9. S.J. Brodsky, Springer Tracts in Modern Physics 100 
(1982) 81. 

10. A.V. Efremov and A.V. Radyushkin, Theor. Math. Phys. 
42 (1980) 97. 

11. G. Bertsch, S.J. Brodsky, A.S. Goldhaber, and J. Gunion, 
Phys. Rev. Lett. 47 (1981) 297. 

12. V.L. Chernyak and A.R. Zhitnitsky, Phys. Rep. 112 (1984) 
173. 

13. E791 Collaboration, E.M. Aitala et al, Phys. Rev. Lett. 
86 (2001) 4768. 

14. R. Jakob and P. Kroll, Phys. Lett. B315 (1993) 463. 

15. H.C. Pauli, Nucl. Phys. B (Proc Supp) 90 (2000) 154. 

16. H.C. Pauli and A. Mukherjee, Int. Jour, of Mod. Phys. 
A16 (2001) 4351. 

17. H.C. Pauli, Nucl. Phys. A705 (2002) 73. 

18. G. Sterman and P. Stoler, Ann. Rev. Nuc. Part. Sci. 43 
(1997) 193. 

19. S.J. Brodsky and G.R. Farrar, Phys. Rev. Lett. 31 (1973) 
1153. 

20. L.L. Frankfurt, G.A. Miller, and M. Strikman, Phys. Lett. 
B304 (1993) 1. 



